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Spontaneous  formation  of  nanometer-scale  composition-modulated  structures  is  a common 
phenomenon  for  III— V and  II— VI  semiconductor  alloys  [ ].  These  structures  are  mostly 
formed  in  open  systems  during  the  epitaxial  growth.  It  is  typical  for  semiconductor  epitaxial 
growth  that  bulk  diffusivity  of  atoms  is  negligibly  small  compared  to  the  surface  diffusivity. 
For  an  alloy  growth,  composition  modulation  in  every  atomic  layer  is  being  formed  at 
the  surface  during  the  growth  of  this  layer  and  persists  when  this  layer  is  buried  by  the 
subsequent  layers.  Modulations  of  composition  from  the  entire  thickness  of  the  epitaxial 
film  affect,  via  long-range  strain  field,  the  surface  migration  of  atoms  [ , , ]. 

For  alloys  grown  on  a (00  l)-substrate,  most  of  observed  composition-modulated  struc- 
tures are  2D  ones  modulated  in  both  [100]  and  [010]  directions.  The  possible  directions 
of  modulation  can  be  explained  by  the  linear  stability  analysis  of  the  homogeneous  alloy 
growth.  If  the  elastic  anisotropy  is  the  dominant  effect  compared  to  the  surface  diffusion 
anisotropy  [ ],  the  instability  occurs  for  the  first  time  in  elastically  soft  directions  [100] 
and  [010].  Then,  however,  the  linear  stability  analysis  does  not  address  the  question  of 
the  selection  between  ID  and  2D  structures.  To  describe  the  final  modulated  structure,  a 
non-linear  theory  is  needed. 

In  our  earlier  papers  [ ] we  considered  the  effect  of  non-linear  coupling  between  compo- 
sitional and  morphological  instabilities  on  steady-state  structures  formed  during  epitaxial 
growth.  It  was  shown  that  at  sufficiently  low  growth  velocities  2D  steady-state  structures 
form,  modulated  in  both  [ 1 00]  and  [010]  directions.  This  effect  is  similar  to  that  for  strained 
islands  in  lattice-mismatched  systems  where  2D  structures  such  as  an  array  of  pyramids 
provide  more  efficient  elastic  relaxation  than  ID  structures  such  as  an  array  of  prisms  [ ]. 
However,  for  lattice-matched  systems  this  effect  is  important  only  for  sufficiently  large  pe- 
riods of  the  resulting  structures,  d > 700  A.  This  is  beyond  typical  scales  where  modulated 
structures  could  serve  as  quantum  dots. 

Among  other  non-linear  effects  which  may  lead  to  the  formation  of  2D  structures  the 
effect  of  species-dependent  atomic  mobility  [ ] is  the  most  plausible.  The  Arrhenius- 
type  behavior  of  the  diffusion  coefficient,  D(T)  — Doexp(—Ea/T),  together  with  the 
dependence  of  the  activation  energy  on  composition  Ea  — Ea(c)  results  in  very  strong 
dependence  of  diffusion  coefficient  on  composition.  In  the  present  paper  we  consider  the 
effect  of  this  dependence  on  steady-state  structures.  This  effect  is  purely  kinetic  and  has 
no  analogs  in  thermodynamics. 

We  consider  the  growth  of  an  alloy  Ai_cBfC  by  molecular  beam  epitaxy  on  an  atom- 
ically rough  surface.  The  alloy  is  lattice-matched  on  average  to  the  (OOl)-substrate. 
The  growth  proceeds  via  deposition  of  atoms  on  the  surface,  surface  migration  of  atoms 
in  a stress-  and  composition-dependent  chemical  potential,  and  incorporation  of  atoms 
into  the  growing  crystal,  desorption  being  neglected.  Let  the  composition  be  equal  to 
c(r)  — c + <p( r),  the  average  composition  being  c — 1/2,  and,  for  simplicity,  the  surface 
profile  be  taken  to  be  flat  h(x,  y)  = vt,  where  v is  the  average  growth  velocity  controlled 
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Fig.  1.  Steady-state  diagram  containing  the  regions  of  homogeneous  growth,  of  the  growth  of  a 
ID  structure,  and  of  the  growth  of  a 2D  structure. 


by  the  deposition  flux.  The  kinetic  equation  describes  the  evolution  of  composition  fluctu- 
ations 0 (r)  at  the  advancing  surface 
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Here  F = Fchem  + Fgrad  + Telast  is  the  total  Helmholtz  free  energy,  where  Fchem  ~ 
/ /chem  (c)dV  is  a chemical  free  energy  of  the  alloy  depending  locally  on  the  alloy  compo- 
sition, Fgrad  is  a gradient  energy,  ~ k f (V0)2t/ V,  and  felast  is  the  elastic  energy  induced 
by  composition  fluctuations.  The  composition-dependent  diffusion  coefficient  D(0)  is 
related  to  the  substitutional  diffusion  of  alloy  components  on  the  surface  which  is  taken 
into  account  up  to  the  second  order  in  0,  D(<p)  = Do  + D\(p  + Z>202  and  a is  the  lattice 
parameter. 

The  linear  stability  analysis  [ , ] showed  that,  for  a given  growth  velocity  v,  the 
instability  occurs  in  a certain  temperature  interval.  The  positive  formation  enthalpy  of  the 
alloy  and  the  elastic  interaction  energy  between  the  adatoms  and  the  “buried”  composition 
modulation  in  the  entire  thickness  of  the  epitaxial  film  are  the  two  driving  forces  of  the 
decomposition.  At  high  temperatures,  the  entropy  contribution  to  the  alloy  free  energy 
hinders  the  driving  forces  to  decomposition  and  stabilizes  the  growth  of  a homogeneous 
alloy.  At  low  temperatures,  the  surface  diffusion  is  “freezed”  out  and  does  not  provide 
formation  of  a modulated  structure.  The  solid  line  in  Fig.  1 depicts  the  boundary  of  the 
instability  region  found  in  [ ].  The  homogeneous  growth  is  unstable  at  v < vc(T)  — 
[4Do(Tc  — T)]/((i2k),  where  /’,  is  the  critical  temperature  in  the  slow  deposition  limit  [ ]. 

In  the  present  paper  we  seek  the  steady-state  solution  to  kinetic  equation  (1)  in  the 
weak  segregation  regime  close  to  the  boundary  of  the  instability  region,  i.  e.  at  v/vc  — 
I — i]2,  where  r\  <5C  1-  In  this  regime  there  is  only  one  unstable  mode  for  each  elas- 
tically soft  direction,  and  composition  fluctuations  can  be  written  in  the  form  </>(x,y)  — 
2t]  [(rc  — T)  /a3]  *^2  [0.r  cos  (kox)  + 0V  cos  (koy)]  + stable  modes,  where 
ko  — y/ 2{Tc  — r)/(fl3/c).  The  reaction  of  stable  modes,  such  as  those  with  k = (2ko,  0) 
and  k = (ko,  ko),  on  unstable  modes  can  be  taken  into  account  in  the  adiabatic  approxima- 
tion. 

To  illustrate  this  reaction,  we  consider  a composition  fluctuation  with  the  wave  vector 
k = (ko,  0)  (Fig.  2(a)).  This  fluctuation  is  unstable  and  its  amplification  implies  the  surface 
flux  of  atoms  B from  A-rich  domains  to  B-rich  domains,  Fig.  2(b).  Assume  that  A-rich 
domains  have  larger  atomic  mobility  than  B-rich  domains,  D(<j>  < 0)  > D(<j>  > 0), 
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Fig.  2.  Stabilizing  effect  due  to  the  dependence  of  surface  atomic  mobility  on  composition  in  ID 
system,  (a):  Unstable  composition  fluctuation,  (b):  Surface  flux  forming  composition  fluctuations 
of  Fig.  2(a).  (c):  Surface  atomic  mobility,  (d):  Surface  flux  of  Fig.  2(b)  in  case  of  varied  atomic 
mobility  (c).  (e):  Second  harmonic  of  composition  fluctuations  due  to  modulated  surface  flux  (d). 
(f):  Surface  flux  induced  by  the  second  harmonic  (e).  (g):  The  first  harmonic  of  the  surface  flux  in 
case  of  varied  atomic  mobility  (c). 


Fig.  2(c).  In  this  case  the  flux  of  atoms  B from  A-rich  domains  towards  domain  boundaries 
is  larger  than  the  flux  from  domain  boundaries  to  B-rich  domains  Fig.  2(d).  This  results 
in  the  accumulation  of  atoms  B at  domain  boundaries,  i.e.  in  formation  of  a composition 
fluctuation  with  the  wave  vector  k = (2Ao,  0),  Fig.  2(e).  Such  a fluctuation  is  energetically 
not  favorable  and  decays  via  surface  flux  out  of  domain  boundaries,  Fig.  2(f).  Due  to 
composition  dependence  of  atomic  mobility  this  flux  is  also  deformed  in  a way  that  most 
of  atoms  B move  towards  A-rich  domains.  Fig.  2(g).  As  a result,  there  appears  non-linear 
stabilizing  effect.  This  is  the  action  of  the  composition  fluctuation  along  x-direction  on 
itself  by  means  of  the  second  harmonic.  There  is  a similar  action  of  the  composition 
fluctuation  in  the  x-  direction  on  fluctuations  in  the  y-direction  and  vise  versa  due  to  the 
appearance  of  composition  fluctuations  with  wave  vectors  like  k = (ko,  ko). 

The  resulting  kinetic  equations  in  the  weak  segregation  regime  have  the  following  form, 
very  similar  to  that  of  Ref.  [ ]: 


= h2  [Vu  — avi>\  ~ tfrV'.vV'2]  + 0(r/4) 

^ = h2  \ify  - aL'I'y  - «7'^yV'2]  + 0(lj4)  , (2) 

where 

_ D2  (Tc-T)  D]  IOC Tc-T) 

Dq  A4a3  Dj  9A4a3 
£ D2  2 (Tc-T)  D\  2 f Tc-T)s\Ra 3 

ClT  = 6 X 5 1 7 X t-  X 

Do  A4a3  D\  A4a3  e^Ra3  + (Tc  - T) 


(3) 
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Here  eo  is  the  lattice  mismatch  between  pure  binary  constituents  of  the  alloy,  AC  and  BC,  R 
is  a certain  combination  of  elastic  moduli,  for  GaAs  R — 0.2eV/A3,  and  A4  is  related  to  the 
fourth-order  contribution  to  the  chemical  energy,  (A4  — 6 (94/chem/34c)  at  c — 1/2).  If 
a 1 < a-r,  the  steady  state  structure  is  ID  and,  if ul  > ay,  the  steady  state  structure  is  2D.  It 
follows  from  Eqs.  (3)  that  at  temperatures  close  to  T,  the  final  structure  is  one-dimensional 
and  at  lower  temperatures  (Tc  — T)  > Ra^,  the  final  structure  is  2D. 

This  transition  can  be  explained  as  follows.  The  species  dependent  atomic  mobility  re- 
sults in  “effective  barriers”  hindering  surface  atoms  from  migration  across  domain  bound- 
aries. These  barriers  are  unavoidable  for  a surface  flux  in  case  of  1 D structure  (Fig.  2),  while 
in  case  of  2D  structure  there  is  a possibility  for  a surface  flux  to  go  around  these  barriers, 
e.g.  along  <110  > -directions.  Therefore,  the  effect  of  species-dependent  atomic  mobility 
favors  2D  structures.  The  dependence  of  atomic  mobility  on  composition  enhances  as  the 
growth  temperature  decreases.  As  a result,  below  a certain  temperature  the  steady-state 
structure  is  2D. 

This  effect  is  in  a way  complementary  to  the  effect  of  surface  corrugation.  While  the 
effect  of  surface  corrugation  is  important  in  case  of  a large  mismatch  eo  between  two  binary 
constituents  of  the  alloy  and  at  a large  spatial  scale,  the  effect  of  species-dependent  atomic 
mobility  is  stronger  for  moderate  eo  and  is  essential  at  small  scales.  E.g.  for  eo  = 0.08  the 
final  structure  is  2D  at  periods  50  A and  smaller,  and  for  eo  = 0.02  the  final  structure  is  2D 
at  d < 200  A. 

To  conclude,  the  steady  state  phase  diagram  of  Fig.  1 applies  to  any  alloy.  We  show  a 
possibility  to  form  2D  structures  with  arbitrary  small  wavelength  and  explain  commonly 
observed  2D  structures  modulated  in  [100]  and  [010]  directions  (see,  e.  g.  [ ]).  Kinetic 
phase  transition  between  the  growth  of  1 D structure  and  the  growth  of  2D  structure  provides 
a possibility  to  tune  the  structure  between  ID  and  2D  ones  at  arbitrary  scale  by  varying 
temperature  and/or  growth  velocity. 
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